Mon. Not. R. Astron. Soc. 000, 000-000 (2012) Printed 4 January 2013 (MN lATrfX style file v2.2) 



Gravitational lensing effects on sub-millimetre galaxy counts 
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ABSTRACT 

We study the effects on the number counts of sub-millimetre galaxies due to grav- 
itational lensing. We explore the effects on the magnification cross section due to 
halo density profiles, ellipticity and cosmological parameter (the power-spectrum nor- 
malisation er 8 ). We show that the ellipticity does not strongly affect the magnification 
cross section in gravitational lensing while the halo radial profiles do. Since the bary- 
onic cooling effect is stronger in galaxies than clusters, galactic haloes are more 
concentrated. In light of this, a new scenario of two halo population model is explored 
where galaxies are modeled as a singular isothermal sphere profile and clusters as 
a Navarro, Frenk and White (NFW) profile. We find the transition mass between the 
two has modest effects on the lensing probability. The cosmological parameter a 8 al- 
ters the abundance of haloes and therefore affects our results. Compared with other 
methods, our model is simpler and more realistic. The conclusions of previous works 
is confirm that gravitational lensing is a natural explanation for the number count ex- 
cess at the bright end. 



Key words: cosmology - gravitational lensing - galaxies: clusters: general - sub- 
millimetre: galaxies. 



1 INTRODUCTION 

Sub-millimetre galaxies (SMGs) at high redshift appear to be 
the counterparts of the most luminous star-forming galaxies 
in the local Universe. The high luminosity of these galaxies 
is presumed to be the result of large amounts of star for- 
matio n, 100 - 1000 h' 1 Mq yr _1 and warm dust iBIai n et al.l 
2002). The millimetre and sub-millimetre wavelengths sur- 
veys have provided an important complement to the optical 
and radio searches for distant galaxies. In recent surveys, a 
significant population of hi gh luminosity, high redshift galax- 
ies ha ve been discovered iCoppin et al.ll2006l:[Nearello et al.l 
120071) . The redshift distribution of the SMGs has a nar- 
row distribution with a probable median r edshift of 2 — 3 
(Chapman et al.l2005l:lAretxaaa et al.l2007b . Surveys at sub- 
millimetre wavelengths show that the SMG population has 
a sharp falloff at the bright luminosity end of their luminos- 
ity function. Gravitational lensing by intervening galaxy clus- 
ters and g roups mod i fies the observed number counts sig- 
nificantly 'jBlainlh 996t ILima et al.lfloTobt Uain & L~imall201 it 
IHezaveh & Holderll201 1|) . Recently it has become possible 
to identify the lensed SMGs, e.g. from the grou nd with the 
South Pole Telesco pe (SPT) IVieira et al J 120101) and from 
space with Herschel iGonzal ez-Nuevo et ai1l2012l) . 

* E-mail: xer@nao.cas.cn 



Gravitational lensing probes the cosmology, the mass 
distribution in the univers e, and prov ides a way to study the 
high redshift objects (see iTreulfeoiOl for a review). In the sta- 
tistical study of cosmological gravitational lensing, several as- 
pects have been studied: th e lensing probab i lity of separa- 
tion o f multiple images (e.g. Keeton & Madail liooilll-i et all 
12007b : t he number of giant arcs f ormed by gravitational lens- 
ing (e .g. lBartelmann et al .11 19981: Hi et al.ll2005l:lHoresh efaTI 
1201 1|) and the modifi ed luminosity function (LF) of back- 
ground sources (e.g. ILima et alj 1201 Obi: IHezaveh & Holderl 
1201 iHWvithe et al.ll201 lUWardlow et al.ll201 2l) . 

In this paper, we will focus on the LF of high redshift 
SMGs. There are a relatively large number of SMGs, and 
their redshift is likely high (z > 2), both of which help 
to c reate an excellent so urce population for lensing stud- 
ies (Chapman et al. 2005). To predict the observed source 
counts, the intrinsic LF needs to be known. Furthermore, 
in order to determine the probability (total cross section) 
of gravitational lensing, we need to know several proper- 
ties of lens haloes, e.g. their abundance and internal struc- 
ture, which are affected by the cosmological parameters 
of the universe. More specifically, the cross section due to 
gravitational lensing may be affected by ha lo ellipticity (e.g. 
iRusin & TeomarkteOOlllHuterer et al.l2005l), the radial profile 
of the lens halo iLi & Ostrikerl2002l:IOauri & Keetonl20 04) as 
well as the size of background galaxies iPerrotta et aL !2002: 
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iHezaveh & Holder|[20TlT> . In lLima etall 1201 Oal) , haloes are 
modelled with NFW profiles, while ellipticity is added to their 
lensin g potential to incre ase the lensing probability. More re- 
cently iLapi et alj J201 2h . considered a composite model for 
galactic sized halo, where dark matter is mod elled as NFW 
and the stellar component by a Sersic profile JSersidh 963b . 
It is in fact close to the is othermal prof ile. They also use the 
public code GLAFIC from lOaunl j2010l) to study the effect of 
halo ellipticity on the lensing cross section and find that the 
ellipticity only weakly affects the cross section of the isother- 
mal halo lens. Our analytical results confirm their finding. 

Axisymmetric lensing models offer simplicity in the study 
of lens statistics. Several spherical lenses have been stud- 
ied, such as the Singular Isothermal Sphere (SIS) and 
the Navarro-Frenk-W hite profile (e.g. iLi & OstrikeH l2002t 
Oguri & K eetonll2004l) . Due to baryonic cooling, clusters and 
galaxies will have different mass profiles. Therefore, a com- 
bination of two population halo mass profiles will be stud- 
ied in this work. The statistical effect of lensing will be 
affected by the cooling mass scale M coo i (s ee, e.g., the 
study on multipl e images and image splitting bv lLi & O striker 
l2002blChenll2004) . A cooling mass scale of ~ 10 13 ft" 1 Mq 
has been suggested (se e, e.g., iPorciani & Madaul [20001: 
iKochanek & White! 12001 ). On large scales the mass func- 
tion of dark haloes relates to the primordial fluctuations of the 
universe, which can be characterised by the power-spectrum 
normalisation parameter <t 8 . Thus, the total cross section will 
be sensitive to this parameter as well. Moreover, previous 
studies have shown that different aspects also affect the lens- 
i ng efficien cy, like the substructures withi n the dark haloes 
(Oaur j|2006l), central massive blac k holes iMao & WittJl2012t 
I Li et al. 112013) and external shear jHuterer et al.ll2005l) . Most 
of them probably do not play significant roles, and so they will 
not be considered in this paper. 

We will study the lensing efficiency, i.e. the magnifica- 
tion cross section in this paper. We start from a single lens 
halo for different halo profiles, and present the lensing cross 
section dependence on halo properties. We show that the 
halo ellipticity does not affect the cross section significantly, 
while the halo density profile does. As we mentioned before, 
a combination of halo profiles is more physical and we find 
such a scenario can produce a higher lensing probability than 
a single universal NFW profile. We present our calculation 
and employ it to study the luminosity function in Section 3; 
we further discuss our results in Section 4. The cosmology 
that we adopt in this paper is a ACDM model with parameters 
based on the results of the Wilkinson Microw ave Anisotropy 
Probe seven year data iKom atsu et al.ll201 1|> : fi A = 0.734, 
fin, = 0.266, fi b = 0.0449, n = 0.963, a Hubble constant 
H = lOO/i kms -1 Mpc -1 and h = 0.71. We allow the cr a 
parameter to vary, but use the WMAP7 value a$ = 0.8 if not 
mentioned. 



2 BASIC FORMALISM 

The fundamentals of gravitat ional lensing can be found in 
Bartelmann & Schneider 1 2001 ). For its elegance and brevity, 
we shall use the complex notation. The thin-lens approxi- 
mation is adopted, implying that the lensing mass distribu- 
tion can be projected onto the lens plane perpendicular to 
the line-of-sight. We introduce angular coordinates with re- 



spect to the line-of-sight. The lensing convergence, that is the 
dimensionless projected surface-mass density, can be writ- 
ten as 

c 2 D. 



k(0) = £(0)/£ cr , where £ c 



(1) 



4ttG D d D ds 

is the critical surface mass density depending on the angular- 
diameter distances D e , D d and D ds from the observer to the 
source, the observer to the lens, and the lens to the source, 
respectively. £(0) is the projected surface-mass density of 
the lens. All lensing quantities can be derived from the effec- 
tive lensing potential ij>, 



ip(0) = - [ d 2 6'K(0')\n\0-0'\ 

7T J K 2 



(2) 



To the lowest order, image distortions caused by gravitational 
lensing are described by the complex shear 

1 



7=2 (^i^ - <9 2 ^) + idid2i> ■ 

The magnification for a point source is given by 

1 

A 4 = 
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2.1 Lensing properties of different dark matter halo 
profiles 

Having laid out a general formalism, we will present the basic 
lensing properties of two different dark matter halo profiles 
that will be used in this paper (see below). They are Singu- 
lar Isothermal Sphe re (SIS) and Navarro-Frenk-White (NFW, 
iNavarro et al.lR 997) profiles. 

The dimensionless surface mass density and shear for 
an SIS halo are 



20 
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(5) 



where (j> is the position angle around the lens, 6 = sj6\ + 6\ 
is the angular separation, and 6> E is the Einstein angular ra 
dius, which is calculated by 



(6) 



where <j v is the one-dimensional velocity dispersion, and c 
is the speed of light. In the rest of the paper, c is used as 
the concentration parameter. The l ensing prope rties of the 
NFW profile has been calculated by Bartelm anrc jl996l) . The 
convergence is analytically given by 



2k, 



/(*) 



(7) 



where x = 6D d /r a (or x = 6/6 s ,6 s = r s /D d ) is the dimen- 
sionless radius, and the function f(x) is defined as 

arcsecha; 



/(*) 



1 - 



(as < i); 



(x = l); 



(x>l). 



(8) 



Vx 2 - 1 

The physical properties of the halo are contained in the pa- 
rameter K a = /9 C ritA c r s /£ C r, where A c is the dimension- 
less characteristic density and p„it is the critical density. 
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Figure 1 . Cross section cr(fi > 2) vs. halo ellipticity for the SIE and NIE 
halo profiles. The solid line is the analytical result for the SIE halo given 
by Eq. 1 1 1 The open circles and dashed line are the cross section using 
cell counts in the magnification map, where the magnification is derived 
by Eqs. l9land lA5l for SIE and NIE model respectively. The open squares 
represent the numerical results for the NIE halo which magnification is de- 
rived numerically (Appendix B). (9e = 5.5 arcsec is used here for all cases 
(which corresponds to the halo mass of M 20 o = 10 14 h- 1 Mq). For the 
NIE halo, a core radius 9 C = 2 arcsec is used. 



The halo mass is defined as M 20 o = 800/3p C rit7rr2oo, where 
r2oo = r s c and c is the concentration parameter (see the ap- 
pendix in Navarro et al. 1997). M ore lensing properties of t he 
NFW halo profile can be found in Wright & Brainerd I2000). 



2.2 Lensing cross sections of dark matter haloes 

We will calculate the cross section of different halo profiles, 
and present their dependence on different parameters, e.g. 
halo ellipticity and concentration etc. Some previous studies 
have shown that the magnification cross section incre ases 
with halo ellipticity dramatically (e.g. lLimaet al.l l2010a). We, 
however, find that this is mainly because the halo mass is 
also changed. We start with an elliptical surface density for 
which the mass within an ellipse is identical with that of the 
axis-symmetric one. In general, the corresponding m agnifi- 
cation can be calculated numerically iSchramml ll990|) . 

We first perform a series of tests for the cross section 
of the halo ellipticity. The halo ellipticity is characterised by 
e = 1 - q, where q = 9 a /9 b is the axis ratio, 9 a and b are 
the major and minor axes, respectiv ely. The an gular sepa- 
ration 9 will be replaced by 6 = \J q9\ + 6f/q (uppercase 
6 is used for the elliptical coordinate) in Eq.[5]to calculate 
Kj and 7 for a Singular Is othermal Ellipsoid (SIE) halo (see 
Keeton & Koch aneklll998l for more details on the lensing 
properties of an SIE halo). It is a radial coordinate and will 
be constant on elliptical contours (thus x = &D d /r s in the 
NFW model). The total mass within a given 6 is invariant 
with the ellipticity e. Following Eq.[A5] the magnification of an 
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Figure 2. Differential lensing cross section as a function of magnification 
fj, for the SIE halo. $e = 5.5 arcsec is used as in Fig.[TJ The solid and 
dashed lines are the analytical and numerical results, which almost overlap. 
The dot-dashed line shows the expected asymptotic behaviour. 



SIE halo has an analytical expression 
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For a given magnification jLt m in, the cross section o-(^) is the 
area inside which the magnification of a source is equal to or 
larger than ^ min on the source plane. For the SIE halo, the 
cross section can be given analytically 

1 1 



+ 



(10) 



>-l) 2 (m + 1) 2 
Notice that at large magnification, the cross section fol- 
l ows the predicted as ymptotic power-law l/fj 2 for n > 1 
iSchneideret al.lll992l) . 

The mass of the SIE halo is related to 9 E - For the SIE 
model the mass is related to the rotation vel ocity, which ca n 
be calculated using M 200 = V?/(WGH(z)) iMo et al.ll 998). 
The relation between the rotation velocity and velocity disper- 
sion is complicated and may be different for different types of 
galaxies. We use the approximate relation between the rota- 
tion velocity an d the velocity dispersion a v = Vc/V%, which 
is suggested bv lChad I2010I) . The cross section has an ex- 
pression of mass and distance 



(T(fJ,) 



16ty j 



2/3 



(GHM 2 . 



x4/3 



D s 



+ 



( M +l) 



(11) 



In our definition of the ellipticity, the mass and the cross sec- 
tion do not change with ellipticity. We can see that only the 
halo mass and the redshifts of lens and source affect the 
cross section. In the appendix, one can find more analyti- 
cal results for the magnification and cross section of SIE and 
Non-singular Isothermal Ellipsoid (NIE) profiles. 

To check the accuracy of Eq.Qo] we create a magnifica- 
tion map using Eq.|9]and sum up the pixels following Eg.lA6l 
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A halo with mass of M200 = 10 14 h' 1 M Q is used, which 
gives 6> E = 5.5 arcsec for z d = 0.5 and z B = 2.0. The res- 
olution of the magnification grid is 6e/500. We increase the 
lens halo ellipticity to see its effect to the cross section. For 
simplicity, we use a constant halo ellipticity as a function of 
radius. The numerical results are shown by the open circles 
in Fig. [Hand the red solid line shows the prediction of Eg.fTol 
We can see that the results agree well with each other. It 
encourages us to explore whether such an independence is 
also valid for other kinds of mass distribution models. Fig. [2] 
presents the probability of cross section as a function of mag- 
nification for the SIE halo {da/d/x). The solid line represents 
the theoretical prediction, while the dashed line is obtained 
by numerical calculation. One can see that the numerical re- 
sult closely follows the theoretical prediction. As expected, at 
small n the probability density distribution rapidly decreases; 
at large /x, it becomes close to the asymptotic relation oc n~ 3 . 

By adding a core of 8 C = 2 arcsec into the SIE model, 
we extend our test to the NIE model. In Appendix A, we show 
an analytical expression of the magnification (Eq. |A5). By 
summing up the pixels in the magnification map as before, 
we show the "theoretical" cross section as a dashed line in 
Fig. [T] We call it the theoretical cross section because the 
magnification is calculated analytically. We also show the re- 
sult where the magnification is calculated numerically. One 
can see that the cross section of NIE halo increases with 
the ellipticity by a few percent for the parameters considered 
here. 

For most elliptical mass distribution models, the 
magnification c an not be derived analytically , e.g., the 
Einasto model (Retana-Montenea roet al.N20la . the Hern- 
quist mod el iBaes & Deionghel feoo^l) and the NFW model. 
Schramm (1990) proposed a way to calculate the lensing 
properties for any kin d of elliptical m ass model. This algo- 
rithm is summarised in lKeetonll200ll) and revised for our def- 
inition of elliptical coordinate in Appendix B. Following this ap- 
proach, we recalculate the magnification numerically for the 
NIE halo and plot the results in Fig. [T] as the open squares. 
The good agreement between the dashed line and the open 
squares guarantees the accuracy of the numerically calcu- 
lated magnification and confirms the weak dependence be- 
tween the cross section and the ellipticity. 

For an elliptical NFW (eNFW) halo, the lensing proper- 
ties, e.g. shear, magnification and the cross section, can be 
calculated numerically as above. We perform similar numer- 
ical tests for the cross section a(fi) using eNFW halo pro- 
files. The same mass M 20 o = 10 14 h' 1 M© and redshifts 
(z d = 0.5, z s = 2.0) are used. 

In Fig. [3] we show the cross section <r(/x > 2) varia- 
tion with halo ellipticity. Two different concentration parame- 
ters c = 5,7 are used for NFW halo profiles. We find that 
there is a weak dependence on the halo ellipticity of a (less 
than 5 percent). Fig.|4]shows the cross section variation with 
the concentration parameter c for the NFW profiles. As ex- 
pected, with a higher concentration the lensing efficiency is 
higher. Fig. [5] shows the cross section dependence on 
the magnification fx. At high magnification, there is a signif- 
icant difference in cross section between two concentration 
parameters for the NFW profile. 

There are a number of complications we have ignored 
(see also the discussion). The numerical tests we performed 
only take into account the projected shape, and assumed the 
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Figure 3. Cumulative cross section cr(fi > 2) vs. halo ellipticity for the 
NFW halo profiles: The solid (dashed) line is the result for the NFW halo 
profile with concentration c = 5 (c = 7). The same condition, halo mass 
M 2 oo = 10 14 h' 1 Mq and redshift z d = 0.5, z s = 2.0 are used for all 
the figures in this section. 
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Figure 4. Lensing cross section vs. concentration for the NFW profiles. 
The solid line is the normalised cross section of fj, > 2 for NFW profiles. 
The dotted line is for the cross section with fi > 4. The same halo mass 
M200 = 10 14 h' 1 Mq is used. 



ellipticity is constant as a function of radius. In reality, this 
may not be true. In particular, irregular shape lens haloes 
may have different cross sections and probability of gener- 
ating multiple images especially for the massive haloes not 
long after the merging event. Other properties, such as the 
substructures may also affect the cross section. The finite 
source size may smooth the magnification and lowers the 
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Figure 5. Lensing cross section vs. fi for the NFW profile. The solid and 
dashed lines represent the results for the NFW profile with c = 5, 7. Again 
a halo mass M200 = 10 14 h~ 1 Mq is used. 



lensing efficiency (see the end of Section 3). Thus, in reality, 
the lensing halo may have slightly different cross sections, 
which need to be further studied using numerical simulation, 
although the required spatial and mass resolutions may be 
challenging. As the effect of ellipticity on the lensing cross 
section seems to be relatively small, in this paper, for sim- 
plicity we will adopt the spherical lens model in the following 
calculations. 



2.3 Probability function of magnification 

The probability distribution P(n) can be estimated either by 
ray-tracing simulations (e.g. Hilbert et al. 2007), or by semi- 



analytical methods (e.g. lLima et al.ll201 b), which integrates 
all the halo contributions along the line of sight from us to 
the source redshift. We will adopt the latter approach, which 
allows us to easily test the dependence on parameters, e.g. 

(j 8 and M coo i. 

Different halo density profiles will be employed for lens 
haloes, i.e. SIS and NFW. The mass-concentration relation 
for th e NFW halo profile has been studied by several authors 
(e.g. lBullock et al.ll200lllZhao et al.ll2009l) . In this paper, we 
adopt the simple model 



c(M,z) = — — I — 



M 

M, 



(12) 



where A/* is calculated by cr m (M») = S c . cr^(M) is the vari- 
ance of the linear density field (Eq.[T5), and S c is the linearly 
extrapolated density contrast threshold at redshift z in spher- 
ical collapse, here we use S c (z = 0) = 1.686. 

The halo mass function jPress & Schechter] 119741) . 
which determines the number of haloes given a mass at each 
redshift n(M, z), can be written as 



AM M 



(13) 



where p is the comoving mean matter density of the U niverse 
and v = 6 c /a m (M). In the lSheth & Tormenl 119991) formal- 
ism, we have, 



A\i —av 

7T 



[1 + (av 2 )~ p ] exp [-au 2 /2] 



(14) 



Here o 2 m (M) is the variance of the linear density field in a 
top hat of radius r that encloses M = 4Tvr 3 p m /3 at the back- 
ground density 



(2ir)3 



\W(kr)\ 2 P L (k), 



where P L {k) is the linear power spectrum 

P L {k,z) oc k n D 2 (z)T 2 (k), 



(15) 



(16) 



and W(kr) is the Fourier transform of the top hat window 
function. The fitti ng formula of the linea r transfer function in- 
cluding baryons (Eisenst ein & Hulll998l) is used here a nd we 
use the WMAP 7-year cosmology iKo matsu et al.ll201 it) . Al- 
ternative models for the shape of n(M , z) are available in 
the literature (e.g. Matarrese et al. 2000), we will not consider 
these since Eq. {14} provides a good description of the mass 
function in numerical simulations. We however use another 
form of Eq. (T4l which is easier to implement in numerical 
integration 



It has a fitting formula of 



M An 
p dlriCT -1 



(17) 



f(a) = 0.315exp(-|lna" 1 + 0.611"^) (18) 

in the range -1.2 < lna -1 < 1.05 jJenkins et aT]l200lh . 

The lensing cross section <r(> fi, z d , M) for a single halo 
will be calculated as a function of mass. The halo mass is 
used as A/ 20 o. The cross section of the SIS halo can be cal- 
culated analytically (Eq.QoJ and that of the NFW halo will be 
performed numerically. The sum of all cross sections in the 
Universe can be written as (Lima et al.ll2010al> 



CTtot (A*) = 47T 



dM AM M >fr**M). 

(19) 

where D A is the comoving angular diameter distance, and 
we make a simplification that the source SMGs are dis- 
tributed at a fixed redshift. We place th e sub-millimetre 
source galaxies mos tly at redshift z s = 3.0 I Chapma n et al.l 
2005: L^n et al.ll200"7l) . but will later allow it to vary between 
2-4 (see Section 3). The lens haloes are distributed between 
the sources and us (0.001 < z d < 1.5). 

The total cross section is integrated using Eq.[T9] then 
the cumulative probability that a source at z s is magnified by 
a factor greater than is then 



p{n) = 



Qtot(AQ 

47T 



(20) 



The probability density can be obtained by p(/x) = 
—AP(p)/Au. The magnification probability function will be 
affected by several aspects, e.g. the lens halo profile, halo 
mass function, and the lens and the source redshifts etc. We 
will study several factors. First of all, the cross section of a 
single halo <r(> a,z d ,M) can be calculated using different 
models. We perform the calculation for three profiles: 1) the 
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SIS halo profile, 2) the NFW halo profile, and 3) a two popu- 
lation combination of SIS and NFW profiles where the transi- 
tion occurs at the cooling mass scale M coo i between galaxies 
and clusters. 

In the left panel of Fig. [6] one can see that the lens halo 
profile will affect the probability function p{y). All the results 
follow p(fi) oc fi~ 3 as expected when fj, > 1. At small /i, all 
the profiles generate similar lensing probability and all the 
curves drop rapidly with increasing ^. The SIS model will 
generate a larger probability for large u; when u approaches 
to 1, the SIS model has a slightly smaller probability than the 
NFW profile. The result using SIS+NFW however is close to 
that of the NFW profile. The reason is that the main difference 
between SIS and NFW is for large and the main contribu- 
tion to large « is from massive haloes, i.e. the NFW haloes 
here. 

As expected, the transition mass A/ coo i also affects the 
probability function. For a larger M coo i, there will be more 
lenses modeled as SIS, as a result one will obtain a higher 
probability at large u (middle panel in Fig. [§). In particu- 
lar, we find that the probability is about 10% lower if we 
use A'/cooi = 10 12 h' 1 M Q at large n than the case with 
Mcooi = 10 13 h' 1 Mq ; at small fi there is a few percent dif- 
ference. In the right panel of Fig.[H] one can see that a larger 
as increases the lensing probability due to a larger number of 
massive structures. The difference exists for most \i, although 
it is more significant for large n ( ~ 25%) than for smaller fi 
(~ 10%). 

In addition, we use different source redshifts. The same 
lens halo redshifts are used (0.001 < z d < 1.5), since at 
high redshift, the number density of lens halo will dramati- 
cally decrease. In Fig.0 we can see that the lensing prob- 
ability increases with the source redshift; at large redshift it 
increases more slowly. At very high redshifts (e.g., z s > 5.0), 
the luminous source number density is low, thus high redshift 
sources will not strongly affect our result. Moreover, different 
halo profiles have different source redshift dependence. The 
probability of NFW profiles increase faster than that of SIS 
profile. 



3 NUMBER COUNTS OF SUB-MILLIMETRE GALAXIES 

The intrinsic number density distribution of a population of 
gal axies can be fitted by empirical or semi-anal ytical mod- 
els i Bauah ~et"aTll2005l) . In this paper we adopt a lSchechteil 
j 19761) form for the intrinsic luminosity function: 



dn 
dS 



n 



-s/s* 



(21) 



where n*, S* , and a are free parameters. Lensing by inter- 
vening haloes changes the intrinsic dn/dS to its observed 
counterpart. In addition to making sources appear brighter, 
gravitational lensing also dilutes the source number den- 
sit y by magnifying th e observed solid angle. As discussed 
in lJain & Lima! 1201 1|) , the observed number counts for the 
whole population is 



dn 1 f 1 dp dn' , , . ,, , 
dS = TJt) J /L? df? dS" ^ = ' ' 



(22) 



"mm = 1-1 to ^ max = 30, and use probability (l-P )<5(/i = 1) 
for /i < u min , where P is the cumulative probability in image 
plane from u min to /i max , and 5(u = 1) is the Dirac function. 
We discuss how the choice of « min affects the results at the 
end of this section. 

First, we compare our prediction with the results given in 
the right panel of Fig. 2 in lLima et ail d201 Obi) . We used the 
rescaled intrinsic Schechter function 



dN' _ la _ s , 
dS> ~ ' 



(23) 



The probability at small magnification is difficult to calculate 
using the halo model. We perform the integral (Eq[22) from 



where N' = n/n* and 5" = S/S*. Different n* and S* are 
used f or different wavelengths (for more detail see I Li ma et al.l 
2010b). We adopt a combination of spherical NFW and SIS 
halo profiles for lens and show our prediction in Fig. [8] The 
rescaled luminosity functions before and after lensing are 
shown by the lines and the points are the rescaled data us- 
ing different wavelengths and different surveys. One can see 
that our r esult can match the data as well as the model pre- 
diction of lLima efaTI 1201 Obi) . From the view of methodology, 
our model is simpler and more meaningful. Their model has 
to adopt a high ellipticity of halo (e.g. 0.4) to reach the re- 
quired lensing efficiency. We also show the effect of the tran- 
sition mass Mcooi- We allow M coo] to vary from 10 12 h' 1 M 
to 1O 14 /i _1 M . The predicted uncertainty is shown by the 
shaded region. The model with a larger M coo i will predict a 
higher number count at the bright end due to more effective 
SIS halo lenses. 

We then apply our method to calculate the sub- 
millimetre galaxy counts of HerMES observation. We use 
n* = 5 x 10 3 /deg 2 , s* = 10 mJy a nd a = -1.0 ob - 
tained by fitting the low flux data points iGlenn et aTll2010l) . 
In Figs. [9] and [TO] the lines show the luminosity function be- 
fore and after lensing. Here we mainly consider a source red- 
shift z s = 3.0. The SIS profile will generate more lensed 
images tha n the NFW profile. The po i nts sh ow the Her- 
MES data I HerME S~C~ollaboration et alj 1201 2d to compare 
with our predictions. Similar as Fig. [6] different models are 
compared. Lensing does not significantly affect the galaxy 
number counts at low luminosity, while at the bright end, lens- 
ing dramatically enhances the number counts. But all the pre- 
dictions with our model (sources at redshift z s = 3.0, a 8 = 0.8 
and M coo i = 10 13 h~ x M Q ) have lower number counts than 
that observed by HerMES. Wardlo wet al.l I2012T) also inves- 
tigated the lensing effect on the HerMES SMG data and find 
that this discrepancy is mainly due to the contamination from 
the late-type sprials. Although our lensing efficiency has been 
enhanced by the SIS halo, the discrepancy is still significant. 
The number of lensed S MGs in which thes e sprial galaxies 
have been subtracted I Wardl owet al.ll2012l) is shown as the 
solid points in the figure. Our results confirms that gravita- 
tional lensing is a natural way to explain the behaviour of 
SMG luminosity function at the bright end. 

In Fig. [9] the shaded region shows the uncertainty due 
to source redshift between 2.0 and 4.0. All of our predictions 
with a single source redshift are larger than the lensed data 
(solid points). In reality, the source will have a redshift distri- 
bution. Different source redshift distributions are suggested 
due to different selection criteria. A significant number of low 
redshift sources a re also found in the Herschel survey (e.g. 
Case y et aTll2012l) . Furthermore, the Schechter luminosity 
function parameters such as S* or a may be different as well. 
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Figure 6. Lensing probability density p(^t) along the line of sight up to z B = 3.0 from all intervening haloes. Left: Results using different halo profiles of 
SIS, SIS+NFW and NFW are shown as the red solid, green dashed, and blue sh ort dashed lines res pectively. M coo i = 10 13 h~ 1 Mq is used. The black 
dot-dashed line represents the combined probability of log-normal approximation l lHilbert et al.l201 lb and SIS+NFW. Middle: Results for the SIS+NFW halo 
profile using different A/ cool 's: 10 13 ri _1 Mq (green dashed line), 10 12 h~ x Mq (blue short dashed line). Right: Results for SIS+NFW halo profile using 
different erg's: erg = 0.7 (red solid line), erg = 0.8 (green dashed line), erg = 0.9 (blue short dashed line). In all three panels, the green dashed line represents 
the result using the same parameters (SIS+NFW halo, M CDD i = 10 13 /i" 1 Mq and erg = 0.8) for better comparisons. 



So a more realistic prediction will depend on these unknown 
parameters, nevertheless the agreement is encouraging. 

In Fig. [To] we show the effects due to cosmological pa- 
rameter cr 8 . As one expects, a large <t 8 will increase the halo 
number density for all masses, which will increase the lensed 
galaxy number counts at both bright and faint ends. From the 
shaded regions in Figs. [8] [9] and [To] one can see that the 
uncertainty of source redshift affects our result the most. 

In addition, we perform our calculation using different 
^ min (1.05 and 1.15). The effects on the number counts are 
small, 5% a t the b right end and 7% at the faint end. In 
Hilbe rt et all 1201 a log-normal approximation of lensing 
probability is found using large volu me numerical simu lations. 
We use the fitting formula given by Hilbe rt et ai1i201 1|) to cal- 
culate the probability for fj, < 1.1. We perform the integral 
from fimin = 0.01 to /i max = 30 and find that the change 
to the number counts at the bright end is also small (blue 
line in Fig.fTOl The upper limits /Lt ma * is uncertain as well. 
iPerrotta et alj J2002fa consider the SMG population and esti- 
mate /imax ~ 10 - 30. However, a magnification of u ~ 45 is 
reported in the cluster Abell 221 8 iKneib et al.ll2004l) . Our re- 
sults show that it causes a decrease of 50% at the bright end 
{S > 0.1 Jy) if we adopt = 20 but keep the faint end un- 
changed. This effects will be blurred if the finite source size 
of background galaxies is taken into account and therefore 
needs a more detailed study. 
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Figure 7. Lensing probability density p(fj, = 2) from intervening haloes 
(0.001 < < 1.5) for different source redshift z s . The solid, dashed and 
dotted lines represent the probabilities for the SIS, NFW and SIS+NFW 
profiles respectively. 



4 SUMMARY AND DISCUSSION 

In this paper we have studied the lensing effects of galaxy 
and cluster haloes, and how lensing will change the lumi- 
nosity function of high redshift sub-millimetre galaxies. The 
lensing properties for individual haloes and haloes as a pop- 
ulation have been presented. In particular, we find that the 
halo ellipticity does not affect the lensing efficiency signifi- 
cantly with our definition of ellipse coordinate. On the other 



hand, as expected, halo mass profiles significantly affect the 
lensing cross section. Not surprisingly, the NFW profile has 
lower lensing efficiency than the SIS profile. We argued that a 
combination of two population of halo profiles is more realis- 
tic: we use the SIS model for mass less than A/ coo i and NFW 
for mass greater than M cool . The SIS profile is favoured for 
galactic sized haloes, due to the baryonic cooling effect. It is 
also close to the composite profile of an NFW da r k halo and 
a Sers ic stellar component iGavazzi et al.l i2007l) ; lLapi et al.l 
(2012). Our prediction with simple assumptions (two pop- 
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Figure 8. The scaled intrinsic and lensed galaxy number counts dN' /dS': 
the solid red line is the intrinsic Schechter luminosity function. The black 
line with grey shaded region represents the lensed galaxy number count due 
to the combination of halo profiles for sources at redshift z s = 3.0, and 
the shaded region displays the range of lensing predictions due to different 
transition masse s M coo i. Th e open circles are data points from the right 
panel of Fig. 2 in lLima et al] tolObh . 



ulation, circular sy mmetric halo profi le) can match the ob- 
servation given in lLima et ail 1201 Obi) quite well. Compared 
with previous works, our model is simpler and more realistic. 
For th e SMG number counts of (HerMES Collaboration et al. 
|2012|) , the excess on the bight end can not be explained by 
lensing alone. The contamination from local late-type sprial 
galaxies is non-negligible. On the other hand, the cleaned 
lensed SMGs sample can be easily fit by lensing. We vary 
M coo i between 10 12 and 10 13 h' 1 M , and find it only slightly 
affects the lensing probability. Moreover, our results also 
show the dependence of the lensing probability on the cos- 
mological parameter, <t 8 . This is especially important for large 
magnifications, since in the halo model, the number density 
of massive haloes is sensitive to cr 8 . 

In a simple test of total cross section as a function of halo 
mass, we find that most contribution to the lensing cross sec- 
tion comes from dark matter haloes with mass between 10 12 
and 10 15 ft _1 M Q , i.e. massive galaxies or galaxy groups. 
The upper mass limit is due to the rare number of massive 
haloes. It also explains why most lenses are at low redshift, 
since massive haloes have not yet formed at earlier times. 
The lower limit is due to the small cross sections of individ- 
ual low-mass lens haloes, because they generate little cross 
sections. 

Our model is simplistic, a number of improvements can 
be made. (1 ) We did not take into account the additional mat- 
ter distribution between the source and observer, i.e. sev- 
eral small haloes may contribute to a modest lensing mag- 
nification. The contribution t o the cross section appears to 
be small, e.g. a few percent jHilbert et aTll2007l) . (2) We as- 
sumed the ellipticity is constant as a function o f radius. For 
relaxed haloes, this may be a fair approximation | Jing & Sutol 
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Figure 9. Intrinsic and lensed galaxy number counts: the red thick solid 
line is the intrinsic Schechter luminosity function describing galaxy number 
counts. The other lines represent the lensed galaxy number counts due to 
different lens halo profiles: SIS (cerulean dot-dashed line), NFW (green 
dashed line) and SIS+NFW (dark solid line), assuming that all the source 
galaxies are at redshift z s = 3.0. The shaded region displays the range of 
prediction due to different source redshifts (from z s = 2.0 to z s = 4.0). 
The poin ts are cumulative 500 fim number counts for H erMES blank-field 
catalogs jOliver et al|2010l:lHerMES Collaboration et alj|2012h and P(D) 
analysis iGlenn et alj2010h . 



2002). Ho wever, baryons may make the central parts more 
spherical iSpringel et al. 2004). Even more significantly, for 
merging galaxies or clusters, the cross section may be en- 
hanced. (3) We ignored the fi nite source size, which m ay re- 
duce the magnification effect iHezaveh & Holder||201 Us - 
ing a smaller maximum magnification (Perrotta et al. 2002), 
the predicted number count becomes smaller at the bright 
end. However, the overall effect to the lensing probability due 
to finite source size needs a more detailed study. The lim- 
ited resolution of Herschel may also confuse multiple images 
as a single one, and thus change the magnification bias. (4) 
The multiple images due to substructures and external shear 
may increase the lensing efficiency. A more detailed study 
with numerical simulations is desirable to address some of 
these issues. 
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APPENDIX A: CROSS SECTION OF SIE AND NIE 
HALOES 



From lKeeton & Kochanekl |1998|) . the surface density for the 
NIE profile is 

bi 



n(x,y) = 



Its magnification reads 



2^q 2 (x 2 + s 2 ) + ;. 



H = 1 r + 



bis 



(A1) 



(A2) 



ip + s) 2 + (1 - q 2 )x 2 } ' 

where V 2 = <7 2 (s 2 + x 2 + y 2 ). Now we define an elliptical 
coordinate, = \/q0f + 0%/q, and write the surface density 
as 

0E 



k(0) = 



2ve 2 + e 2 ' 



(A3) 



These definitions keep the total mass within an ellipse invari- 
ant with ellipticity, q for given Qe,0 c and 6. We can rewrite 
Eq.|A3las 



(A4) 



Comparing Eqs. |Al] and |A4] the magnification in our defini- 
tion is 

// : i(-)..r; - 1 



* *[g(* + 6» c /g) 2 + (l-q 2 )0 2 ]' 



(A5) 



where * = a/9 2 + 9 2 . The cross section for a given magnifi- 
cation threshold is defined as 



f(Mmin) = 



d 2 (3 = 



/J|>/Jmin 



pr d 2 6> (A6) 



and it can always be calculated by numerical integration. The 
SIE model has C = and has an ellipsoidal magnification 
distribution. The corresponding cross section can be derived 
analytically (see Eqs. |9] and [To}. 



|2001|) . Scaling coordinates as x = ^/q9\ and y = s/qd 2 , we 
can rewrite the ellipsoidal distribution as 

k = «(9), where 6 2 = qOl + Ol/q = x 2 + y 2 /q 2 . (B1) 

The lensing properties for a surface density distribution with 
elliptical symmetry can be written as a set of one-dimensional 



(B2) 

(B3) 
(B4) 
(B5) 
(B6) 
(B7) 



integrals, 




ip(x,y) 


= §/(*,v) 




= qxj (x,y) 




= qyJi(x,y) 


1pxm{x,y) 


= qx 2 K (x,y) + qj () {x,y) 


i>vv(x,y) 


= qy 2 K 2 {x,y) + qJ 1 (x,y) 


ipx V (x,y) 


= qxyKi(x,y). 


Here 
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(B8) 
(B9) 
(B10) 



are one-dimensional integrals, where V-v is the circular de- 
flection angle, q is the axis ratio of the lens and «' are the first 
order derivatives of the convergence, i.e. «'(£) = d/«(^)/d^. 
The convergence k is written as a function of the ellipse co- 
ordinate £(u) given by 



x 2 + 



1 - (1 - q 2 )u 



(B11) 



Therefore once we know the behavior of k'(t-) and Vv(r) of 
a circular symmetric «;(r), the lensing properties of the ellipti- 
cal mass distribution «(6) can always be obtained through 
these integrals. The above equations are slightly different 
from those in lKeetonll200ll) . Considering the scaling relation 
between coordinates (a;, y) and (9i, 6 2 ), we simply have 



tpBiBi (01,02 

ipe 1 e 2 {0i,02 
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(B12) 
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APPENDIX B: NUMERICAL METHOD FOR ELLIPTICAL 
HALO LENSING PROPERTIES 



The lensing properties of an elliptical halo can be cal culated 
numerically given an arbitrary surface density profile iKeetonl 



